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Atom interferometer phase in the presence of proof mass 

B. Dubetsky 
hdubetskyQgmail. com 
(Dated: May 6, 2015) 

This is presented the justification of the expression for the atom interferometer phase in the pres¬ 
ence of proof mass nsed in http://arxiv.org/abs/1407.7287. Quantnm corrections to that expression 
are also derived. The corrections allow one to calculate numerically atom interferometer phase with 
accuracy Ippm or better. 


I. MAIN RELATIONS 

For times between the Raman pulses the atomic density matrix evolves as 

ih'p = [H, p ], (1) 


where 


H = 


2Ma 


U (x, t) 


( 2 ) 


is Hamiltonian, p, Ma and U are, respectively, atomic momentum, mass and gravity potential. In the Wigner repre¬ 
sentation 


p(x,p,t) = — J dsp(^x + ^s,x-^s,t^ exp{-ip-s/h) 


one 


finds [I|, Sec. Ill 


( 3 ) 


dt + jfds - dsU (x, t)dir+Ql p (x,p, t) = 0, 


Q = - [ih) 

When the size of the gravity potential 


1 


U [ X —ihdp, t \ — U [ X — -ihdp, t 


1 


-I- dsU (x, t) dp. 


(4a) 

(4b) 


L^U[x)/\dsU (x)| 


( 5 ) 


is sufficiently large so that 


ApL 


< I 


( 6 ) 


where Ap is a size of the momentum dependence of the density matrix, one expand over the dp term and gets 
approximately 

Q ~ -^Xikiix,t)dp,dp^dp^ (7a) 

Xikiix,t) = -doo.doo^d^^U {x,t) (7b) 


An implicit summation convention of Eq. © will be used in all subsequent equations. Repeated indices and symbols 
appearing on the right-hand-side (rhs) of an equation are to be summed over, unless they also appear on the left-hand 
side (Ihs) of that equation. For the relative weight of Q—term one finds 


Q 

dsU (x, t) dp 


h? 

24 [LApf 


« 1 , 


so that one can consider Q—term as a small perturbation and accept that 


( 8 ) 


P {x,p, t) = Po (x,p, t) + Pq (x,p, t ), 


(9) 
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where {x, p, t) is unperturbed density matrix obeying equation 


dt + - dsU (f, t) dp Pq {x,p, t) =0 


and Pq {x,p,t) is a perturbation evolving as 


dt + - dgU (f, t) dp Pq {x,p, t) = -Qp^ {x,p, t) 


(10) 


( 11 ) 


We assume that density matrix is known at some preceding time t', i.e. Eqs. (nnmn) are subject to initial conditions 


PQ{x,p,t') = p{x,p,t'). 
PQ{x,p,t') = 0. 

Solution of the homogeneous Eq. m is given by [I|, Sec III 


Po 1) = Po (-R 1) - P 0 > 1') > 


(12a) 

(12b) 


(13) 


where |i? {x,p, <i, <2) , P {x,p^ ii, i2)| are atomic classical position and momentum at time ti subject to initial condi¬ 
tions {x^p} at time t 2 - Atomic classical trajectory, evidently, obeys the multiplication law. 


(^R{x,p,t',t") ,P {x,p,t',t'') =1^1 {x,p,t,t"). 


(14) 


Consider now solution of the Eq. (HU. Since Q—term is initially equal 0, solution of the inhomogeneous Eq. (HU 
contains only the stimulated part consisting of contributions produced at different times t' < t" < t. Since operator 
inside the curly brackets in Eq. HU is a full time derivative, 

^ =dt + -^ds - dsU (f, t) dp 
dt -‘.Via 

in the time interval + dP'] one produces contributions to the Q—term 

dpg {x,p, t" -P dt") = -Qpq {x,p, t") dt". (15) 

Later on this contribution evolves freely [i.e. satisfies Eq. HU subject to the initial condition HU at time t" + dt"\ 
and at time t one finds for this contribution 


dpQ {x,p,t) = -^dt" 


Xtki dstds^dtf^po 




Total solution is a sum of all these contributions 


Pq {x,p,t) = — dt" 






Substituting here solution (TT^ for Pq one finds 


Pq ix,p, dt" Xiki (I t") {P (C, t', t”) , P t") , t') 


^ ^ —-j ^ ^ {x,p,t" ,t) 


(16) 


(17) 


(18) 


II. AI PHASE 


TT TT 

Consider now an atom cloud launched at time to , and interacting with '2 '2 Raman pulses applied 


at times 


T — {to + ti:to + ti + T, to + ti + 2T} , 


(19) 
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where ti is time delay between cloud launch and 1st Raman pulse. T is time delay between pulses. We assume that 
Raman pulses produce coherence between two hyperfine sublevels g and e of the atomic ground state manifold, and 
that initially atomic density matrix m is given by 


Pgg ix,p,to) = f{x,p), 

Pen {x,p,to) = Pgg {x,p,to) = 0. 


(20a) 

(20b) 


One finds, for example in [I|, that after ——pulse applied at time t, density matrix elements jump to the values 


Pee ^ + 0) = ^ Pee ^ “ 0) + Pgg (x,p- Hk, t - o) 
Pgg {x,p, ^ + 0) = ^ Pee {x,P+ hk, t - 0^ + Pgg {x,p, t - 0) 

( 


Peg ix,p,t+0) = - exp 


+ Re < i exp 


— Re < i exp 


hk 


—i{U-x — Si 2 t — 


—i{k-x — 5i2t — 


Peg 1 1 

(21a) 

/ hk ^ 

Pen ix,P+—,t-0 


Pee h?..P+ - 0 -Pggl X,P- ^ 


k ■ X — 6i2t — 

+ ^{Pes(^^.P>^-0) + exp 2i (k ■ X - Si2t - (j)j p^g (x,p, t - 0)| 
and after tt— pulse, the density matrix elements jump to the values 

Pee ix,P, t + 0) = Pgg (x,p- hk, t-Oj , 

Pgg ix,P, t + 0) = Pgg (x,p+ hk, t-0^ , 

Peg (x,p, t + 0) = exp 2i {kx - 5i2t - Pg^ {x,p, t-0). 


hk 


(21b) 

(21c) 

(22a) 

(22b) 

(22c) 


where k is effective wave vector, 5i2 is detuning between fields’ frequency difference and hyperfine transition frequency, 
(j) is phase difference between traveling components of the Raman field. We allow pulses to have different detunings 
and phases (5l2^ (* = 2,3). 

Our purpose is to obtain the atomic density matrix of the excited state Pgg {x, p, ra + 0). We will achieve this by 
applying consequently Eqs. (mini) for density matrix evolution between Raman pulses and before the 1st Raman 
pulse, and Eqs. dm [22]) for density matrix jumps after the pulse. 

Before the first pulse, the density matrix becomes 


After this ——pulse 


Pee {x,P,Ti+ 0 ) 
Pgg {x,p,Ti+ 0 ) 

Peg {x,p,Ti+0) 


Pgg {.X,p,Ti - 0) = / (^R{x,p,to,Ti),P{x,p,to,Ti)'^ 


(23) 


i/ (^R {x,p-hk,to,Ti'^ ,P{x,p- 

i/ [R{x,p,tQ,Tx) ,P {x,p,tQ,Tx)^ , 


-2 exp 


i{% ■ X — P12X1 — I 


0 


J ^ ^ hk 
I I A I a:,p- —,to,Ti 




(24a) 

(24b) 

(24c) 


One uses this density matrix as an initial value at moment ti for the free evolution between 1st and 2nd pulses of 
the unperturbed density matrix 


Po {x,P,Tx +0) = p{x,p,Ti +0) 


(25) 
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We now consider Q—term (fTSl) before the second pulse action. From Eqs. (ITSl I24cl [25l) one obtains 

t2 pT2 


r 

PQeg ix,p,T2-0) = J 


dt 


X 


/ 


exp 

i(k- R (|,7f,Ti,t) - 



' 



R (^R 

\ \ / ^ /j 

(7 -> \ ^ - \ hk ' 

\ 



^ Xikl 

X/ 

,P (^5,7r,Ti,tJ - —,to,Ti 

(^,TT,Tl,t) ,P (^,Tr,Tl,t) - 




P ( R 

< 


\ V 

V J \ > ^ , 

1 


y 


(26) 




In the Eq. (1261) . the derivative on the momentum if is of the order of Ap where Ap is the momentum size of the 
factor in brackets. This factor has two terms, phase-factor and atom distribution /. Initially, for t = ti phase-factor 

is if—independent (because R if, ti, and derivative is of the order of 


where 


^^^Thermal ^ Pq ; 


po = 


is thermal momentum, the Tc is cloud temperature. 


Consider now the phase factor exp 


i{k-R 


at t — Ti 


estimates lets ’’turn off’ the gravity field. Then atom trajectory evidently equal to 


R 


= i- {t-Ti) 


TT 

~a 


(27) 


(28) 


T. For the purpose of the 


(29) 


and the phase factor acquires Doppler phase k - - (t — ti) . It becomes a rapidly oscillating function of momentum 

TT having period of the order of 


Pd ^ Ma/kT 


and phase factor derivative is of magnitude 


O — 1 

^^iDoppler PD • 

For the ^^^Cs at temperature Tq ~ 3pK, k = 1.4743261770* 10^m“^,T = 160ms one finds 

da i Thermal 1 


2* 10“® < 1, 




Doppler 


kvnT 


(30) 


(31) 


(32) 


where vq = po /Ma is thermal velocity. Condition (1321) means that time separation between pulses T is sufficiently large 
to be in the Doppler limiting case, when Doppler phase kv^T 3> 1. In this limit atomic distribution over momentum 
is sufficiently smooth to neglect its derivation, and throughout the text we include only contribution to the Q—term 
arising from the derivative of phase factor. For this reason we did not consider above the Q—term for cloud evolution 
before 1st pulse. As we will show, the atomic levels’ populations (p^g and Pgg) have no phase factor at tg < ^ < '’" 3 ) 
and therefore Q— term arises only from atomic coherence p^g. Using the Doppler limiting case ()32p results in 
AI phase pretty much independent of the atomic momentum and spatial distribution. 
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Calculating derivatives, one obtains 


PQeg {x,P-,r2 “ 0) = J dt 

^u^v^wXikl ^TZi^u j I'^li^ ^TZi ^1 j 




exp 


i{k- R (^|,7r,ri,t) - 


X/ 


^ rIr ,P (^^,7f,Tl,t) - ^ 


P R 


,P (^,7r,ri,t^ - ^,to,Ti 


hk 

~ 

hk 

” 


/ 


(33) 


^ ^ |■{s,p,t,T2) 


In the described approximation, the derivative no longer acts on the term inside the curly brackets. Applying the 
multiplication law M we find that 



7r,Ti,t 


' i=R{S,p,t,T2),Tf=P(x,p,t,T2) 


R 

P 


(x,P,Ti,T2) . 


(34) 


The expression inside the curly brackets of Eq. (1331) becomes t—independent and the Q—term in front of the 2nd 
pulse is then given by 


PQeg (.X,p, T2-0) = |exp i {k ■ I - “ <(>l) / (C, } 


R{x,'P,T\ ,T2) 


^ J ( P(x,p,Ti ,T2) — ^k/2 


pT2 

X dt 

Jti 


kutkeuX^kl Ti, Rv (c, 7?, r 1 , Rnj (i, 7?, n , 


5 r~i ^ 


(35) 


From Eqs. (UlIMlIIl 

Pee {x,P,T2 - 0) = i/ (^R (^^,7f,to,ri) , P TT, tg, Ti) ) 


^ — R{x,P,Ti ,T2),TT — P{x,P,Ti ,T2) — hk 


Pgg {x,P,T2 - 0) = i/ (^R{x,p,to,T2),P{x,P,to,T2)j , 
Poeg ix,p,T 2 “ 0) = -^ |exp i (k 


x/(p (|,7f,to,Tl) ,P(C,7f,to,Tl))} 


^ — R(x,P,Ti ,T2),‘7^—P{x,p,Ti,T2) — hk/2 


(36a) 

(36b) 

(36c) 


From Eqs. (I22|) after the 2nd pulse density matrix jumps to the value 
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Peeix,P,T2+0) = ^f(R(S,p-hk,to,T2^,P(x,p-hk,to,T2^^, 

Pgg{x,P,T2+0) = i/ (|,7r,io,Tl) ,P 


^—R(^x,p-\-hk,Tl,T2^^'^—P(^x,p-\-hk,Tl ,r2^—hk 


Poeg ix,P, T2 + 0) = ^ |exp |i ^k ■ (2x -t) - ‘2S%^T2 + 5^12T1 - 24>2 + 4>^ I 

X/ {R ,-P 


^—R{x,p^t\ ,T2),^—P{x^p,T\,r2)—hk/2 


(37a) 

(37b) 

(37c) 


PQeg{x,P,T2+0) = - | CXp | Z |^fc • - 2(5^2^ T2 + ^^2^ Ti - 2(;i2 + <(>i | 


'-f{R (C,’f,^o,Ti) ,p(|,7r,<o,Ti))| 


^—R{x,p^t\ ,T2) ■,^—P{x^p,T\,r2)—hk/2 


j7T,to,T 1 

kukykwXikl 
dsiRu (^^,Tf,Ti,t^ djf^Rni 


X dt 
J T 1 


^ r“i ^ \{.s,p,t,T2) 


(37d) 


Consider now Q—term produced during free evolution inside the time interval [t 2 , T 3 ]. Each density matrix element 
(1371) produces Q—term. However, since diagonal matrix elements (I37al I37bl) contain no rapidly oscillating in momen¬ 
tum space phase factors and we neglect their Q — terms. Term p7dl) is already linear in Q term and can produce only 
higher order non linear in Q contributions, which we are not including yet. So one has to consider only the Q—term 
produced by coherence (I37a . which we denote as p'qeg- From Eq. (fTSl) one finds 

P'Qeg{x,p,Tz-0) = J { XiM i , CXp |i 2k ■ R (j,n,T2,tj 


'T 2 

- k-R 




- 2(5)2^T2 S\yTi - 2(j)2 + (/fij I 
x/ [r (ji (^R (^i,n,T 2 ,tJ ,P (|,7f,r2,t) ,ri,T 2 ) , 

(■R (c, T2, t) , p {I, 7?, r2, t) , Ti, r2) - Y’ *0, Ti 

[R (c, 7?, 7‘2, t) , P (c, 7f, r2, t) , Ti, T 2 ) , 

{R (l,'^,T 2 ,t) ,P {l,TT,T2,t^ ,Ti,T2^ - ^Rq,Ti 


P 

P 

P 


(38) 

^ ^ = '1 j>(x,p,t.T 3 ) 


where we used the multiplication law da, 

P {r (^^,7r,T2,t) ,P{l,^,T2,t^ ,T1,T2) = P 7?, Tl, . 


(39) 

In Eq. (1381) we differentiate over momentum tt only the rapidly oscillating exponent. After differentiation, we apply 
the multiplication law two more times to the phase factor and distribution /, namely 


(C,7f, r 2 ,t) 




{x,P,T2,T3) , 


(40) 


and therefore 


(P [l7f,T2,t) ,P(e , TT, T2, , Tl, T2^ 

{ ^} (-^(^>^’'^ 2 , 7 - 3 ) ,P(f,p,T2,T3) ,Ti,T2) ^ 




(41) 
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to find out that these terms become t—independent. As a result one gets for Q—term p'q^^ before the 3-rd pulse 
PQea - 0) = jexpjz k ■ (2 R{x,p,T2,T3) - - 2 S^^^T 2 + - 2 (j )2 + 


x/(^ (C,7?,to,Ti) ,.p(|,7r,to,Ti))| 

(^Ru (i,TT,Tl,t^ - 2 Ru (^i,^,T 2 ,tJJ 


^ ^ 


X / ‘fxzfc/ kukyk^ 

J T2 

(^Rv Tl, t) - ‘2Rv (i, TT, '^2, t) ) 

dsi (^Rw (c, TT, Tl, - 2R^ TT, r 2 , ^ I 




(42) 


Pee(.X,P,T3-0) = 

[R (C.7?Po,T2) 

,p(, 

C.tr,to,T2)) 

Pgg{x,p,T3-0) = 


l,p(, 



Applying Eq. (fT^ one hnds that other matrix elements (1571) in front of the 3rd pulse become 

^—R{x,p,T2,T3),7T — P{x,p,T2^T3) — hk 

{i=R{R{x ,P,T2,r3),P{x,p,T2,T3)-\-hk,Ti,T2), 

7T = P(^R{x,p,T2,T3),P{x,p,T2,T3)Phk,Ti ,r2)—^fc} 

P0egix,P,T3 -0) = ||exp |i 2fcP (f, p, T 2 , T 3 ) - fc • | “ 2(5i2^T2 + T i - 202 + </>l } 

Xf (R(t^,to,Ti) ,p (t^,to,ri))]^ ^ ^ , 

V V /V // ^ ^—R{x^p^Ti,r^),'K—P{x,p,Ti,Ts) — hk/2 

PQeg {x,P,T3 - 0) = |exp |z |^2fcP (f, p, T 2 , T 3 ) - k ■ - 2S‘i^ T 2 + - 202 +0l] } 

x/ (a ,P TT, to,Tl)') V ^ 

V \ / \ // ) ^—R{x,p,Ti,T^),7T — P{x,p,Ti,T^) — hk/2 

ku^V^w\ikl ^TTiRu dj^j^Ry 7r,T3,t^ 


(43a) 

(43b) 


(43c) 


pr2 

X dt 






(43d) 


Aggregating different parts of the nondiagonal density matrix elements (I43cl I43d[ HTl) one concludes that in the linear 
in Q—term approximation this term acts only on the atomic coherence phase: 


Peg {x,p, T 3 - 0) Ri (f,p, r3 - 0) Pg^g {x,p, T 3 - 0) + p'g^g (f,p, T 3 - 0) 


^{exp{* 




2kR (f,p, T 2 , Ts) - fc • ^ - 0Q {x,p) - 2^1^72 + - 202 + 0ij I 

\—R{x,p^r\ ,r3),7r—P(x,p,ri,r3)—?).fc/2 


x/ (.R (|’,7f,to,Tl) ,P(^C,7f,to,Tl))| 

'^q[x,p) = 24 


(44a) 


^2 I' ^ 7 * 2 

-—kjkykw \ I dtx'iki {^p') 9siRu dskRv ds,Rw tt, n, 


dtXiki 

' ' 

7 ?, Ti,t^ - 2 Ru 


[Rv[ 

i,TT,Tl,t^ - 2 Ry 1 

|, 7 r,T 2 ,<))j 


(^Rw ( 

|, 7 r,Ti,t^ - 2 Ry, 

(«'"’"“-‘))]}l 

iH"}‘ 


(44b) 




















For the- tt -AI, after the 3rd pulse, one should calculate only the atomic distribution in the excited sublevel. 

One hnds from the Eqs. (I21al I43al I43b[ Hill 


Peg (f,p,T3 +0) = i/ ,P (C,7f,to,T2)) 

+ ^/ (-R (^,^,^ 0 , 1 - 1 ) ,P (^^,7f,to,Tl)) 


— < cos 

2l 

sis). 


hk 


kx — 2kR x,p -—,T 2 ,T 3 + k ■ ^ + (j)Q x,p — 


^—R{x,P,T2,T3),^—P{x,P,T 2 ,T^)—hk 


^ I I R(R(x,p—hk,T2,T3j,P(^X,p — hk,T2,T3^+hk,Ti,T2) 

^ ) y p(^R(^x,p—hk,T2.,r3^,P(^x,p—hk,T2.,r3^Phk,Ti,T2^ — hk 

hk\ 


-Sl'ps + 25^2^T2 - SYpi -(1)3 + 2(^2 - (t^l 
xf (^R (l,7?Po,Ti) ,p(^,7r,to,ri))| 




^—R(^x,p—hk/2,Ti ,^—P(^x,p—hk/2,ri ^T3^—hk/2 


(45) 


This density matrix should be used to calculate any response associated with atoms on the excited state. We are 
using it to get the probability of atom cloud excitation defined as 


w = y dxdppf.^ iS:,p, T 3 + 0 ). 


(46) 


First two terms in Eq. (1451) are responsible for background. Since phase space stays invariant under the atom free 
motion and recoil during interactions with Raman pulses, the background equals to - and 


w = ^{l-w), 


(47) 


where interferometric term is given by 


J dxdp <1 c 


hk 


kx-2kR I x,p- ^,T 2 ,T 3 ) + fc - (()Q I x,p- ^ + 2S\Yt2 - <y) 2 ^ri - ((<3 + 2^2 “ <(>i 


d2). 




^—R(^x^p—hk/2,Ti ,r3'j,7T—P(^x^p—hk/2^Ti ,r3^ — ?ifc/2 


(48) 


Selecting as integration variables atomic initial position and momentum at the moment of atom launching 


{x,p} = |p {i, if, to, Tl) , P TT, to, Tl) 1 

one obtains 


(49) 



(50a) 

[x,]^ = |P ^P (£',/, Ti, to) ,P(*',.^, Tl, to) + ^^/2,r3,Tij , 


P (^R{x\p,Ti,to),P{x',p,Ti,to) + hk/2,T3,Ti^ 

) + fik/ 2 } , 

(50b) 

\d{x,p} /d{x',f}\ = 1, 


(50c) 

Rix,p-^,T 2 ,T 3 \ = P (^P(f',/,Ti,to) ,P(f',_^,ri,to) + hfc/2,T2,ri^ 


(50d) 


After that, replacing {x',p^} —>■ {x,p} , one obtains 
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where the phase of the AI is given by 


(t>{x,p) 

(t>r {X,P) 
(t>Q {X,p) 


(t>r {X,p) + (j)Q {X,]^ 

k 


R (I, T3,Ti^ - 2R TT, T2, Ti) + ^ 
'^Q ),P(I ,7f,T3,Ti^ 


- {i=Rix,p,Ti,to),Tf=P(x,p,Ti,to)+hk/2~^ 
{^=R{x,p,ti ,to),T!=P{x,p,Ti,to)+flk/2^ 


(52a) 

(52b) 

(52c) 


where (j)Q is given by Eq. (I44bp . 

Part (j)^ includes ’’classical” part of the phase (erasing in the limit H ^ 0) and recoil effect during interacion with 
Raman pulses. For the rotating spherical Earth this this part had been calculated in [ij. We calculate here additions 
to (j)j. caused by proof mass field. 

Part (j)Q is originated from quantum correction Q to the density matrix in Wigner representation equation in time 
between pulses. If atom trajectory is smaler than the size of the gravity field ([S]), one can expand gravity field in the 
vicinity of atoms launching point. Holding only gravity-gradient terms in that expansion one gets (pg = 0, because 
tensor dlEl dissapears. For this reason we did not consider Q—term in [l|. Part (pg one has to know for precize 
measurements of the Newtonian gravitational constant Keeping in mind this application we will calculate pg 

below. 


A. Part 


Propagation functions ,P (af,p,t,t')|, i.e. atomic position and momentum at time t subject to initial 

value {x,p^ at moment t', evolve as 


P{x,p,t,t') 


Ma 


R{x,p,t,t') = 

P{x,p,t,t') = Ma^g + 6g R{x,p,t,t') , 


(53a) 

(53b) 


where g is Earth gravity field, Sg{x,t) is proof mass gravity field. We neglect in Eq. (l5^ the gravity-gradient, 
centrifugal and Coriolis forces caused by the rotating Earth. When Sg (x, t) is a perturbation, approximate solutions 
of the Eq. (ITT|) are Q 


R{x,p,t,t') « R^^'> {x,p,t,t') + SR{x,p,t,t'), 
P{x,p,t,t') Rs P^°^ {x,p,t,t')+6P{x,p,t,t'), 

R^°^x,p,t,t') = X + ^ {t-t') + g ^^ , 

(f,p,t,t') = p +Mag{t-t'); 

SR{x,p,t,t') = dt"{t — t")5g R^^^ {x^pR'Pt') ,t" 

SP{x,p,t,t') = dt"5g R^^^ {x,p,t”,t') ,f 

Functions <5i?, (5P| obey multiplication laws 

(x,p,t',t") ,P^°'> (x,p,t',t") 

ds aproximate expression for propagator . 

Using Eqs. (IMIISSI) one consequently finds 


(54a) 

(54b) 

(54c) 

(54d) 

(54e) 

(54f) 


p(o) 


p(o) J (upUU )■ 

(55a) 

SR\,^ 

J ix,P,t,t ). 

(55b) 


hk \ 

’,P,Tl,to) ,P{x,p,Tl,to) + 

Y,r„ri j . 
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RiR{x,p,Ti,to),P{x,p,Ti,to) + ^,Ts,Ti] \ ^ ^ ^ f^k 

V ^ / \ P^ Hx,p,Ti,to)+SPix,p,Ti,to) + —,Ts,Tl 


i?(°) {x,p, Tl,to) + SR {x,p,Tl,to) , 

hk 

T 


+SR f {x,p,Tl,to) ,P^°'> {x,p,Tl,to) + ^,Ts,Ti 


1 


= R^°'> {x,p,Ti,to) + SR {x,p,Ti,to) + 

^ /■■’'I r 1 ftk 

{x,p,Ti,to) + Ma dt'Sg R^°^ {x,p,t',to) ,t' 


'to 


1 2 
{Ts -Tl) + -g{Ts -Tl) 


+ / dt' {ts — t') Sg 


^(0) 1^(0) {x,p,TiRo) + ,tA ,t' 


Since from Eqs. (1551 154cl) . 

R‘'°Hx,P,Ts,h) = R^°^ (^R^°Hx,P,Ti,to) ,P^°Hx,P,TiRo) ,Ts,Ti^ 

= R^°^ {x,p,Ti,to) + {x,P,TiRq) {Ts-Ti) + l^gp-s - Tif 


one rewrites 


^ - hk \ -» hk 

R {Rix,P,Ti,to),P{x,p,Ti,to) + ^ {x,P,Ts,to) + 2^ {Xs - Tl) 

+ (ts-ti) f dt"5g R^°^ {x,p,t',to) ,t' + f dt” {ti - t') Sg R^°'> {x,p,t',to) ,t' 

Jto Jto 


+ / dt' {ts — t') Sg 


^(0) I ^(0) [x,p,Ti,to) ,P'S^'> ix,p,Ti,to) + 


One uses now that 


^(0) / ^(0) (^^^jp^ri,to) ,P^°^ {x,p,Ti,to) + (t' - Tl) 


x,p,TiRo) ,P^°'’ {x,p \ T 1,to) ,t^ T 


= R^°^ {x,p, t', to) + [f - Tl) , 


and therefore 


(56) 


(57) 


(58) 


(59) 


R \^{x,p,Ti,to) ,P{x,p,tiRo) + ^ {.x,p,TsRo) + ■^- (ts -ti) 

+ (ts-ti) [ dt”Sg R^'^^ {x,p,t',to) ,t' + [ dt”{Ti-t')Sg R^°^ {x,p,t',to) ,t' 


'to 


'to 


+ dt'{Ts-t')Sg R^^Ax,]}^', to) ,t' 


+ / dt' (ts - t') < Sg 


hh 

{x,p,t',to) + (t' - Tl),t' 


- Sg {x,p,t',to) ,t' 


(60) 
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The propagator (I54el) can be rewritten as 

SR{x,p,Ts,to) = [ dt' {Ts-t')Sg ,t' 


'to 


dt'{ ts- t')Sg\R^°^x,p,t',to) ,t'] + / dt'( ts- t')Sg\M°^ {x,p,t',to) ,t' 


'to 


dt' {ts -t')6g R^°'> {x,p,t',to) ,t' + / dt'{Ti-t')Sg R^°'> {x,p, t', to), t' 


'to 


+ {ts-ti) dt'Sg R^°'> {x,p,t',to) ,t' 


'to 


(61) 


which coincides with the sum of 3rd, 4th and 5th terms on the right-hand-side (rhs) of Eq. ([^n|l. Finally 

^ ^ hk \ ^ hk 

R lR{x,p,Ti,to),P{x,p,Ti,to) + I (x,p,Ts,to) +^P(x,p,Ts,to) + (t^ - Ti) 


+ / dt' {ts - t') < 5g 


-» hh 

{x,p,t',to) + (t' - Tl),t' 


- Sg ix,p,t',to) ,t' 


(62) 


Substituting this result in the brackets of Eq. (I52b() for the 1st (s = 3) and 2nd (s = 2) terms, one finds that the 
phase of the atom interferometer consist of the terms corresponding to Earth gravity and the proof mass gravity field, 
and quantum correction, 


(jirix.p) = (j)o{x,p)+5(l){x,p)+(j)q{x,p), 

= R^^'^ {x,p,T3,tQ)-2R^°'^ {x,p,T 2 ,to)+R^°'’{x,p,Ti,to) =k^'^, 
5(j)[x,p) = k- 6R{x,p,T3,to) -26R{x,p,T2,to) +SR{x,p,Ti,to) , 

(t>q{x,p) = k-i>g, 


= [ dt' ix3 - t') 

JtI 


Sg 


-2 / (t 2 - t') 6g 


R^°'> {x,p,t',to) + {f - Ti),t' 

hh 

ix,p,t',to) + (t' -Ti),t' 


-6g R‘^°'> {x,p,t',to) ,t' 
- Sg {x,p,t',to) ,t' 


(63a) 

(63b) 

(63c) 

(63d) 


(63e) 


Using Eq. (I54el) and piecinmg together terms corresponding to the integrations in time intervals [to, ti] , [ti, T 2 ], and 
[r 2 ,r 3 ] one gets 0 


S(j) = k- (T 3 U 30 - tiU 2 o + U 21 -U 31 ), 
Uai3 = / dttd5g(^^d+vt + gE^—,t\, 


(64a) 

(64b) 


Consider now quantum correction (I63dl) . vector (I63el) can be rewritten as 


i’q= dt' {t 3 - t') < 5g 


hh 

Rd^'^ {x,p,t',to) + (t' -Ti),t' 


+ / dt' [t' - Ti) 5g 


R^°'> {x,p,t',to) + {t' - Ti),t' 


-Sg R^°^ {x,p,t',to) ,t' 

- Sg {x,p,t',tQ) ,t' 


(65) 


Substituting t' = T 2 + 0 for the 1st term of Eq. (EID, and t' = Ti + 9 for the 2nd term one finds 


= / d9{{T-9) 


Sg(R^°'> (f,p,T 2 + 9, to) + (T + 6») ,t 2 + 6* j - Sg(^R^°^ {x,p,T 2 + 9,to) ,t 2 + 9^ 


+ 9 


Sg(R^°^ (f,p, n + 0 ,^ 0 ) + +6li -Sg(^3°^ (5,P, ri +9,to),Ti + 0^ 


( 66 ) 
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B. Part (fiQ 


Consider now Q—term (I52b|) . Since 

R \ ( R(^Rix,p,Ti,to),Pix,p,Ti,to) + hk/ 2 ,T 3 ,Ti'j , 

P } yp (^R{x,p,Ti,to) ,P{x,p,Ti,to) + hk/ 2 ,T 3 ,Ti'j ,t,T 3 


(^R{x,p,Ti,to) ,P{x,p,Ti,to) + hk/2,t,Ti 

(67) 


this term becomes 

(t^Q{x,p) = -^kuKkw^j dtxM(i,'t)dsiPu(^t^,xi,t^ds^Rv(i,Tf,Ti,t^ds^Ryj(^^,Tr,Ti,t^ 

+ ( dtx'iki il.t) ds^{Ru{l,^,Ti,t \ -2Ru(l,^,T2Ry\ (^Ry - 2Ry 

J T2 

X dsJRyy(?,n,Ti,t)-2Ryy(^,n,T2,t]] } r -) r -) ■ (68) 

When atoms move between pulses under the action of the homogeneous Earth gravity field g and small but inhomo¬ 
geneous perturbation 6 g{x, t) caused by the proof mass tensor IzEl is caused only by the proof mass, 


Xikl ^aXikh 

Xiki = ^^^^^,Sg^{x,t). 


(69a) 

(69b) 


While we calculate the AI phase in the linear in proof mass gravity approximation, it is sufficient to calculate the 
atom trajectory in Eq. (I68p in the absence of the proof mass, when 


{x,p,t,t') = I I {x,p,t,t') 


and therefore 


Using this result and Eqs. (fTOl ISSl) 


dp.Rj {x,p, t, t') = - {t — t') 


M„. 


(70) 


(71) 


(t>Q{x,p) = J dtXijl[^,tj (t-Ti) 


dtXiji t) [t-s - tp 


^—R(^R{x,p,ti ,to),P(x,p,Ti ,to)-\-hk/2,t,Ti^ 


(72) 


Since we are interested in calculating of the Q—term to the 2nd order in recoil momentum hk, we then can neglect 
recoil in the brackets of Eq. (17^ . We then apply the multiplication law (ITil) and obtain 




pQ ix,p) = dtx^ji]^Rix,p,tRo),t^{t-Ti) +J dtXiji]^Rix,p,t,to),t^iT 3 -t) 

Substituting t = ti + 6 for the 1st term of Eq. CS]), and t = T 2 + 0 for the 2nd term one finds 


(73) 


fp P r 

pQ {x,p) = d9 |0 


R{x,p,Ti + e,to),Ti + e +{T-efxiji R{x,p,T2 + 9,to),T2 + o 




(74) 
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